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This paper is the second in a pair treating a new lattice model for nematic media. In addition 
to the familiar isotropic (/) and nematically ordered (N) phases, the phase diagram established 
in the previous paper (Paper I) contains a new, topologically ordered phase (T) occuring at large 
suppression of topological defects and weak nematic interactions. This paper (Paper II) is concerned 
with the experimental signatures of the proposed phase diagram. Specific heat, light scattering and 
magnetic susceptibility near both the N/ T and // T transitions are studied, and critical behavior 
determined. The singular dependences of the Frank constants Ki, K2, K3 and the dielectric tensor 
anisotropy Ae on temperature as T 



Tj^j, are also found. 



PACS: 64.70.Md, 61.30.Jf, 64.60. Cn, 11.15.Ha 
I. INTRODUCTION AND RESULTS 

In the preceding paper [|| (hereafter referred to as "Pa- 
per I"), we showed that the nematic-to- isotropic transi- 
tion need not proceed through a single first order phase 
transition, contrary to long-held belief ||^. In fact, we 
have shown that this disordering can proceed by a pair 
of continuous transitions, with a novel intermediate phase 
possessing non-trivial topological order. The phase dia- 
gram discussed in Paper I (figure 3 there) is expected to 
be generic for nematic materials since the gauge theory 
from which it is derived depends only on the symme- 
try of nematic media. The present paper is devoted to 
the calculation of characteristic signatures of the phase 
transitions into and out of the topologically ordered but 
physically isotropic phase T discovered in Paper I. The 
quantities we study are the specific heat, light scattering, 
and magnetic susceptibility near each transition. We also 
calculate the singular temperature dependences of the 
Frank constants Ki, K2, and K3 as well as the dielectric 
tensor anisotropy Ae. 

An unusual feature of our model is the presence of two 
distinct, physically isotropic phases, T and /, which are 
separated by a continuous phase transition. Unlike the 
familiar transition between liquid and gas (which are also 
both isotropic fluids), there is no latent heat and no crit- 
ical endpoint between T and /; they cannot be smoothly 
interpolated. The qualitative distinction between these 
two phases, both of which are nematically disordered, 
is subtle: in the topologically ordered phase there are 
few configurations with long nematic defects, and the lo- 
cally coarse grained director field can be unambiguously 
converted into a non-singular vector field throughout the 
fluid. The clearest demonstration of this order would be 
the measurement of the energy per unit length of a defect 
that is externally imposed by boundary conditions. Such 



a direct demonstration of topological order, however, is 
likely to be quite difficult. 

Although one cannot use local probes to distinguish 
T from / deep within these phases, the transition be- 
tween them has very characteristic and measurable prop- 
erties. Our prediction is that under certain conditions, 
one should be able to observe critical behavior by all con- 
ventional means (specific heat singularities, light scatter- 
ing, etc.) between two isotropic phases. Observation of 
such critical behavior would support our scenario. Veri- 
fication of the scaling laws derived below would confirm 
the universality classes identified in Paper I. 

A less dramatic but also surprising prediction of Pa- 
per I is the existence of a continuous phase transition 
between the nematically ordered and the topologically 
ordered states. According to conventional wisdom based 
on Landau theory, the development of nematic order from 
any isotropic state should occur as a first order phase 
transition. As we have shown in Paper I, this need not 
be the case if the physically isotropic phase has topolog- 
ical order, since then the transition can be in the univer- 
sality class of the three-dimensional Heisenberg model. 
Detailed predictions based on our model are presented 
below. 

For convenient reference, we summarize our results 
here. Detailed derivations are found in later sections. 
Unless otherwise indicated, reduced temperature relative 
to a particular transition is denoted by t = {T — T^) /T^, 
and the corresponding correlation length is denoted by ^. 
Subscripts + and — denote quantities pertaining above 
{t > 0) and below {t < 0) Tc- "Above" always refers to 
the more disordered phase, so that the isotropic phase is 
above the topologically ordered phase, which in turn is 
above the nematically ordered state. 

1. Specific heat. The singular specific heat near the 
continuous transitions into and out of the topologically 
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ordered phase varies as 



c(r) = A±|i|-". 



(1) 



Near the // T transition, the exponent a takes the (pos- 
itive) three-dimensional Ising model value ||^ a/ = 
0.1085 ± 0.0075, which implies a divergent specific heat. 
The amplitude ratio A^^/A^ is inverted relative to the 
usual Ising spin model, a result of the duality trans- 
formation discussed in part I. Near the T/N transi- 
tion, a takes the three-dimensional Heisenberg value ||] 
an = —0.130 ±0.021, resulting in a cusp. The amplitude 
ratio A^/A^ takes the corresponding universal value. 

2. Polarized light scattering between isotropic 
states. Near the // T transition, the light scattering in- 
tensity Ixx (q, t) for incoming and scattered light with 
parallel polarizations obeys the scaling law 



(li « 1; 
» 1- 



(2) 



Here a/ and vi are the three-dimensional Ising specific 
heat and correlation length exponents a/ = 0.1085 ± 
0.0075, VI = {2- a/)/3 = 0.6305 ± 0.0025, respectively. 
(The Ising correlation length ^ diverges as t~'^' .) 

3. Depolarized light scattering betv^reen 
isotropic states. Near the I/T transition, the light 
scattering intensity I^y (q, t) for incoming and scattered 
beams with perpendicular polarizations shows only a 
weak non-analyticity at q = 0. That is, Ixy{(l,t) can 
be decomposed into an analytic and a singular part 



J-"^(q,t)=g 



9^ < 1; 
9e> 1; 

(3) 



where fi^^ is a universal function. Note that the ana- 
lytic contributions to the scattering dominate the total 
scattering as |q| ^ 0, i ^ 0. 

The exponent in (|) is (1 - ai)/vi = 1.414 ± 0.006, 
which is greater than unity. The singularity is therefore 
merely a divergence in the second derivative d^Ixy/dq^ of 
the scattering intensity. This singularity is not a promis- 
ing signature of the // T transition, since it is so weak. 

4. Depolarized light scattering near the ne- 
matic phase. Near the N/ T transition, the depolarized 
light scattering intensities I^z (q, t) and ly^ (q, t) obey the 
scaling laws 



hz (q, t) ^ q 



2vt 



-7^^(90, 



x,y. 



(4) 



(The incoming polarization z is along the mean director.) 
The very small anomalous dimension of the spin for the 
three-dimensional Heisenberg model is rjn — 0.02 ± 0.01. 
Matching the form (^ onto the known small- wavenumber 
behavior of liz deep in the nematic and isotropic phases 
(for this purpose, Tand / behave identically) implies 



^-2^i/h(1+2>)h)^ 
^-i/ff (l-2))ff ) 



t < 0, < 1; 

» 1; 
t > 0, < 1- 



(5) 



Approaching the transition from the nematic side (t < 
0), this result implies a relationship between the criti- 
cal behavior of the Frank constants and the dielectric 
anisotropy Ae, since 



c(Ae)2 



and 



c(Ae)2 



K2ql + K^q^z 



(6) 



(7) 



with some constant c (we have taken qy — 0). The form 
(^ then implies both that (i) the Frank constants are 
asymptotically equal upon approaching T/vt (due to the 
isotropy of the scaling law) and that (ii) there is a partic- 
ular relation between the exponents governing the van- 
ishing of the Frank constants and of Ae at the N/ T tran- 
sition. 

These are, indeed, our next two conclusions. 

5. Frank constants. The Frank constants approach 
a common value and vanish as 



K, 



T 



'NT 



(8) 



at the N/ T transition. This can be tested by direct mea- 
surements of the Frank constants, either by light scatter- 
ing (see Item (5)) or by Freedericksz instability measure- 
ments (see Item 8). 

The approach of the "splay" Frank constant Ki to 
equality with the "bend" and "twist" constants K2 and 
K3, however, is very slow. Specifically, we expect 



2:3 



1 



2,3 



|ln(|i|)| 



(9) 



Experimentally, this behavior may look very much like 
Ki being proportional to i^2,3 with a constant prefactor 
of order unity. 

6. Dielectric anisotropy near the nematic 
phase. The dielectric anisotropy tensor Ae vanishes like 



Ae 



1^1 



2/3f, 



as T ^ 



NT' 



(10) 



where Ig] = 0.368 ± 0.004 is the Heisenberg order pa- 
rameter exponent. This result for Ae can be tested either 
through light scattering as mentioned above, or by sim- 
ple optical birefringence measurements. Combining with 
our previous results (Items (4) and (5)) for the Frank 
constants and liz gives 



hz{q,t<0)(xq-'\t\ 



2| i|4/3ff-iyff 



(11) 



With the exact relation 2/3 = iy{d — 2 + rj) that is de- 
rived by combining the Fisher and Rushbrooke scaling 
relations |^, we find eq. (||) for -C 1 (and = 3). 
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7. Magnetic response near the nematic phase. 

Near the N/ T transition, the magnetic response of the 
system is given by 

M„ = xiHc, + Ax(t, H) (^hafip - i(5a/5 j (12) 

where a denotes Cartesian components, n is the nematic 
director (always along H in the topologically ordered 
phase), and the isotropic part of the susceptibility xi is 
analytic through the (continuous) N/ T transition. The 
anisotropic part of the susceptibility Ax, on the other 
hand, obeys 



Ax{t,H) 



H 



(13) 



where A and He are non-universal (i.e., system- 
dependent) constants. The exponent is (f>^ = J^2/r/2 = 
0.62 ± 0.01, where is the renormalization group eigen- 
value of the "spin-tensor interaction Q" for the (n = 3) 
Heisenberg model. To second order in the e = 4 — d ex- 
pansion, Ur « 1.77 ± 0.01. The scaling functions f±{x) 
obey 



a; , as a; ^ 0; 
const., as a; ^ 



(14) 



as X 



where^^-1 = 2{1 + riH)iy"y^ = 1.15 ±0.01. There- 
suit (nWj for implies that the non-linear susceptibility 
diverges as the N/ T transition is approached from from 
the isotropic (T) side as 



'dip 



cx for 



H 



«1, 



with 



72 = 2{I3h - ct'x) 



-0.52 ±0.01. 



(15) 



(16) 



Recall that in this topologically ordered, but isotropic 
phase, M is always along the applied field H. From 
( p^ ) and (|l5|), one deduces a non-analytic (though non- 
divergent) temperature dependence of the linear suscep- 
tibility as the N/ T transition is approached from the ne- 
matic side. 



— ^Xi + ^^X = a + bWf'-, 



(17) 



where a and h are constants. Increasing the magnetic 
field so that H ^ Hc\t\'^^, the singular temperature de- 
pendence ( p^ ) of the non-linear susceptibility becomes a 
singular magnetic field dependence 



dm 



cx 



H 



> 1, 



(18) 



which is readily seen to cross over smoothly to the low- 
field {jj^ « 1) result at ~ H,\t\^--. 

8. Preedericksz instabihty. The Freedericksz in- 
stability 1^ provides a means of measuring the ratio of 
the anisotropic susceptibility Ax to the Frank constants 
Ki. When a nematic is confined between parallel plates 
which have been treated to create boundary conditions 
favoring a particular orientation of the director, it will 
assume that orientation even in the presence of an ap- 
plied magnetic field H that favors a different orientation, 
as long as that field is sufficiently small. Usually exper- 
imentalists choose this field to be perpendicular to that 
favored by boundary conditions. When the strength of 
the applied field exceeds a threshold Hp, however, the 
director rotates away from the alignment favored by the 
boundaries to that favored by the field. This is the Freed- 
ericksz instability. The standard result for the threshold 
field is 



H. 




xO(l). 



(19) 



The precise geometry of the experiment determines which 
Frank Ki appears, as well as the 0{\) factor. Inserting 
the known scaling forms 



K,{t,H)^\t\-"fl 



H 



HM 



(20) 



and equation (|T^) for Ax into (|T^), we obtain 

Hp^L-'+'^f^iL/^H), (21) 

where the scaling function f±{x) has the asymptotic 
forms 



as X 
as X ^ 0: 



oo. 



(22) 



for some constants c, c' and c". These relations imply 
that 

r cx 

Hp cx I = ^(-o.99±o.oi 

l cx e 



for t < 0, L > ^; 
for i < 
for t > 0, L > 



(23) 



Unfortunately, the smallness of r]H makes the dependence 
of Hp on the temperature (and on the anomalous length 
scale in the critical regime) extremely weak in the ne- 
matic phase, and hence difficult to discern experimen- 
tally. In practice, one expects to observe an Hp which is 
nearly constant with increasing temperature, until it be- 
gins to drop in the unusual stretched exponential manner 
indicated in (p3|) . The observation of either this stretched 
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exponential or the apparent temperature independence of 
Hp below Tmt would confirm our theory. 

The constancy of Hp below Tmt should be contrasted 
with the predictions of Landau theory in the case of a 
weakly first-order transition. These predictions are 



Ax(T) a go(T) (X 
K,{T)o^Qlo^\t\. 



(24) 



Given these scalings, (|l^ ) would imply that Hp oc \t\^^^, 
instead of (E3h. 



9 Magnetic susceptibility. The elusive // T transi- 
tion can also be detected by magnetic susceptibility mea- 
surements, albeit with more difficulty. We find that the 
non-linear susceptibility has a singularity at the 

transition: 



'dip 



c\t\ 



analytic. 



(25) 



The rest of this paper is devoted to deriving these con- 
clusions. The results quoted above for the specific heat 
require no further discussion, as they follow directly and 
immediately from what has already been shown (in Pa- 
per I) about the partition function, which determines all 
thermodynamic functions. The other results involve var- 
ious correlation functions. 



II. POLARIZED LIGHT SCATTERING NEAR I/T 
TRANSITION 

Consider first light scattering near the // T transition. 
In general, |^ light scattering in an isotropic material is 
given by Ixx{<i) oc (|e(q)p). This scattering is caused by 
local fluctuations in the (isotropic) dielectric constant e 
which in turn are caused by density fluctuations of the 
various components of the material. For simplicity we 
consider here a situation with just one such density p; 
including more does not alter the conclusions. Expand- 
ing the density dependence of e to linear order about 
the average density, we obtain Ixx{<l) cx (|(5yo(q)p). The 
calculation of the polarized light scattering near the // T 
transition thus reduces to determining the fluctuations of 
the (non-ordering) density p. To calculate these correla- 
tions, we must know how density fluctuations enter the 
Hamiltonian. 

Recall that the basic Hamiltonian defining our theory 
(equation of Paper I) is 



S, -Sj+K U.jUjkUkiUu, (26) 
{ijki} 



where the spins Si are three-dimensional unit vectors on 
the sites of a lattice (cubic for convenience), the variables 



neighbor sites, and the second sum in the Hamiltonian 
runs over elementary plaquettes ijkl. 



We showed in Paper I that the spins are irrelevant at 
small J, so they may be neglected when one incorporates 
density-fluctuation effects into an effective renormalized 
lattice gauge Hamiltonian. There are no symmetry re- 
strictions upon the dependence of the Hamiltonian on p. 
For long wavelengths, the dominant dependence is upon 
the value of p itself, and spatial gradients can be ne- 
glected. A Hamiltonian with the appropriate symmetry- 
allowed coupling is 



(27) 



where U abbreviates the renormalized plaquette term 
from eqn. (26) above. (Recall that U measures defect 
density, so that eq. ( p7| ) describes a couphng between 
defect density and fluid density.) In the lattice model, 
the local density variables Pf^ are associated with the 
plaquettes, so we denote this dependence simply by the 
directional index p,. (Each plaquette is then labeled by 
its normal vector.) Kpi is the renormalized defect stiff- 
ness obtained after integrating out the spins as described 
in Paper I, 7 is the coupling between p and disclination 
density, and Sp^^ = p^^{x) - po- 

For the moment, all we need to know about V{p) is 
that it is a local, smooth function through the transition, 
since p is non-ordering, and it has a minimum near (but 
not at!) the equilibrium density po- A duality transfor- 
mation can be applied to the model implied by Hamilto- 
nian (p7|). As explained in Paper I, the duality transfor- 
mation is completely local in the couplings, so the parti- 
tion function becomes 

^ = j cxpI^^{~Vr{Sp^) + J^{Sp^)aa<Ta+f,)^, 



(28) 



with a "renormalized" potential Vr and density- 
dependent spin coupling J given by 



VR{p^) = V{p^)-{Kn + ^Sp^) 

-l/21n(l-e-4^('-^)) 
J{Sp^) = l/21ncoth(ifi^ + 7%). 



(29) 



In this model, fluctuations in p are driven entirely 
by those of the Ising spins that are dual to the gauge 
field. More precisely, they are proportional to the energy- 
energy correlations of the Ising model obeying the scal- 
ing law (||). To see this, expand Vr to second order and 
J{6p) to first order in powers of Sp. The result for the 
Hamiltonian is 

H « Hisingipo) + h^Sp^ - jR^ 5p^ 

+ A/2j2iSpr (30) 



±1 are associated with links (ij) between nearest- with 



h = Vliipo), A = FA'(po), and 77? - J'(po), (31) 
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where primes indicate derivatives. The calculations re- 
quired to find expressions for ((5/9^) and {5p^5pi^i_i) are 
easy in this approximation since they are Gaussian in- 
tegrals over 5p (fluctuations on different sites are even 
decoupled) . Specifically, 



(32) 



Since this is zero by definition, h = ^R{(Ja<^a+ij) ■ Using 
this expression for /i, 

A^iSpfj^Spfj^') = jR{craaa+f^;aa'Cra'+f^') + ASi_,^>, (33) 

where {x;y) = (xy) — {x){y) is the truncated correlation 
function. Aside from the zero-range second term, the 
right hand side of this equation is simply the energy- 
energy correlation function for the three-dimensional 
Ising model, the well-known behavior of which now im- 
plies the result eq. (||) (Item 2). 



III. DEPOLARIZED SCATTERING NEAR N/T 
TRANSITION 



Next, we turn to depolarized light scattering near the 
N/ T transition. The standard form for such scattering 
in an anisotropic medium is |^] 



/y(q) oc (e,j(q)ey(-q)), 



(34) 



where eij(q) is the Fourier transform of the position- 
dependent local dielectric tensor and the summation con- 
vention has been suspended on the right-hand side. 

For small nematic order the dielectric tensor can be 
expanded in the tensor nematic order parameter 



g,,(r) = 5,(r)5,(r)-Jy/3 



Cij — eoSij + gQij, 



(35) 



(36) 



where g is a microscopic parameter, and is therefore non- 
singular at the transition. 
The correlation function 



a,{r) = (e,,(r)e,,(0)) 



(37) 



for i ^ j — z can therefore be expressed in terms of 
Heisenberg spins of our model as 



CUr) oc (S,(r)S,(r)S,(0)S,(0)). 



(38) 



In Paper I we have shown that near the N/ T transi- 
tion the system is described by an effective Heisenberg 
model, so that we can apply a renormalization group 
transformation to relate Ciz{r) to the same function in 
the Heisenberg model with different parameters: 



^9±{r/0- (39) 



The second equality follows upon making the choice b = 
r/a, where a is a microscopic length and ^ = CoN"'^" 
is the three-dimensional Heisenberg correlation length, 
noting the isotropy of the Heisenberg model, and making 
the definition 



g{x) = C,,(a,xi/''-(eo/a)i/''-)a2(^-2+^-). 



(40) 



Fourier transformation (and the substitution d = 3) leads 
directly to eq. (Q) (Item 4), with 



fNTD 

J± 



(PR 



i?2(l+r,) ' 



(41) 



To complete the description of the critical scattering 
near the N/ T transition requires matching the scaling 
form (^ onto known long-wavelength behavior in the N 
and T" phases to obtain the asymptotic forms of f±. 

Consider first the topologically ordered ( T) phase {t > 
0). No Goldstone modes are present, since there is no 
broken symmetry. The fluctuations of the dielectric ten- 
sor (hence light scattering) must therefore be finite as 
q| ^ at fixed temperature. Consistency with the scal- 



ing law demands f+{x) oc x 



as a; ^ 0. The g 



scattering therefore diverges upon approaching the tran- 
sition from the isotropic side: 

/,^(g-^0,i) oce^-2rm ^^^u(i-2^nH)^ qi<^l. (42) 

The requirement of finite scattering at fixed q in the crit- 
ical regime (q^ ^1) implies that /+ and /_ tend to the 
same constant as a; — ^ cxo; this means 



/^.(9) « 9" for ge>i. 



(43) 



On the nematic side of the transition {t < 0), scat- 
tering is dominated by Goldstone modes — fluctuations 
in the orientation of the director. The resulting long- 
wavelength behaviors of eqs. (^, 0) are standard results 
Q. In these equations, the anisotropy of the dielectric 
tensor, Ae, is given by 



Ae 



(44) 



the mean director n lies along z, and q • y = 0. If the 
scaling results are to coincide with eqs. (|[ 0) in the 



q —>■ limit, then f-{x) oc x' 
result, 



-(l+2»)) 



as X 



0. As a 



Irziq) - ^-(i+2riH)q-2 (45) 

OC as ^ and t 0". (46) 

This is the result announced in equation (|^) of Item 4. 
Implications for the critical behavior of the Frank con- 
stants are also carried by eq. ( |4^ ) combined with eqs. (|^, 
Specifically, 



^M!oc|<r-(i+2''-), i = 1,2,3. 
Ki 



(47) 
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Clearly all three Frank constants have the same critical 
behavior. 

We now verify both eq. (|4^) and the critical behavior 
{Ki ^ Itl"") of the Frank constants by an independent 
derivation. If we take an expectation value and express 



(48) 



Qij in terms of spins (see eq. (pSD) we find 
(Qy) oc |<|2/'(n,n, -%/3). 



Comparing with the definition of Ae, we deduce Ae ex 

Some information about the Frank constants can be 
obtained from generalized Josephson relations Q . Imag- 
ine imposing a twist on the boundary conditions for a 
nematic confined to an L x L x L volume, and compare 
to one with periodic boundary conditions. The (Frank) 
free energy difference is AF{L) = {K2/2)L9^ , where 6 is 
the twist angle. From scaling, one expects AF — g{L/^) 
near a continuous transition. Combining these two equa- 
tions. 



\t\' 



(49) 



Very similar arguments that use boundary conditions 
which impose a bend or splay distortion imply Ki , Kt, cx 
It]" as well. Gathering these results and the scahng of 
Ae we find 



(50) 



Application of the exact scaling relation 2/3 = i/((i — 2 + 
77) — v{l + rf) (the last equality following from d = i) 
recovers the result eq. (E^. 



IV. CORRECTIONS TO SCALING FOR FRANK 
CONSTANTS 

Let us next consider the corrections to scaling for the 
Frank elastic constants. In the context of an |S|'' soft- 
spin theory Q , the starting point is a Landau free energy 
expressed in terms of a vector- valued field S with fluctu- 
ating magnitude and direction. The term lowest order in 
gradients is uniquely of the form /(|Sp)(VS)^, as long 
as the invariance under global rotations of S is unbroken 
(compare the discussion in Paper I). In real nematics, 
this invariance is broken and only the invariance under 
combined identical rotations in internal and real space 
survives. 

The Landau- Ginsburg free energy then contains three 
types of gradient terms: 



^grad - i^l(V • S)2 + i^2|S ■ (V X S)|2 4 

i^alS X (V X S)|2, 



(51) 



where Ki, K2 and are the Frank constants for splay, 
bend, and twist, respectively. In an e = 4 — d expansion, 
the renormalization group eigenvalues of these operators 



at the n = 3 Wilson-Fisher fixed point are A2.3 = —2 -|- 
0(e), as can be seen by power counting (four powers of 
S and two gradients). Since these eigenvalues are large 
and negative, the difference between K2 and Kj, vanishes 
rapidly as the transition is approached. The standard 
form for the correction to scaling is 



(52) 



with A2 and A3 given above. 

The difference between K2 and therefore vanishes 
rapidly, but Ki is another story. The eigenvalue of the 
operator (V • S)^ is § Ai = -eVlOS + 0{€^), so this 
Frank splay constant vanishes very slowly under renor- 
malization. By way of illustration, for e = 1 and Ai = 1, 
Ki still differs by 50% from the value i^^ol^T taken by 
K2 and at the ridiculously tiny reduced temperature 

\t\ « 2-100/'' « 10-43! 

In fact, Ai is so small that for all reasonable val- 
ues of Ai, the renormalization group flows of AKi = 
{Ki — K2,'i) are presumably dominated by a non-linear 
term proportional to {AKif' (which to our knowledge 
has never been calculated). The full recursion relation 
for AKi is then of the form 



dAKi 
d£ 



(53) 



where c is a constant (again, to our knowledge, never 
calculated). Due to the smallness of Ai, eq. (|5^ ) is valid 
outside an extremely narrow range of temperatures in 
which AKi ^ Ai. Assuming c > - i.e. the fixed point 
AKi = is globally, as well as locally, stable ~ and ne- 
glecting the linear term, we conclude that AKi vanishes 
as 1/i for all experimentally observable temperatures. 

If we choose £ such that e^^ = 1, we find £ cx | In \t\\~^. 
This in turn implies that {Ki — K2.3) / K2,3 vanishes like 
I ln(|i|)|^i, as claimed in Item (5). This dependence is so 
weak that in many situations it will simply look as if the 
correction term is independent of t. In this case, we can 
write Ki{t) = Ko\t\''{l + Ai) for all practical purposes. 
That is, Ki will vanish with the same exponent as K2 and 
K3, but with a prefactor of (presumably) order one. This 
means that for q • y = 0, /j,z(q) will remain anisotropic: 



4z(q) = 



(1 + Al)g2+92 



(54) 



The singular character of the anisotropy should be unob- 
servable for all practical purposes. /j,z(q) is still expected 
to become isotropic. Note that the anisotropy of I^z does 
not alter the scaling law Ixz oc \t\''^-^^^^''^ q^'^ , but rather 
introduces a (nonsingular) anisotropic prefactor. 



V. MAGNETIC SUSCEPTIBILITY NEAR N/T 
TRANSITION 

The application of a magnetic field H adds to the 
Hamiltonian a term 
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-AxoE(H-S.)', 



(55) 



where the molecular diamagnetic anisotropy A^g is as- 
sumed to be a smooth function of temperature, in the 
spirit of modern critical phenomenology. This standard 
expression eq. (p5h for the nematic magnetic energy is 
the lowest order J1CI| ] term in H which respects the gauge 
symmetry (Si — > cnSi, Uij (TiUjUij). 

The magnetization M induced by the applied magnetic 
field can be deduced via the standard thermodynamic 
relation 



1 dF 



NdHa 
Now we know that 



(56) 



(57) 



where Qo{T,H) is the magnitude of the nematic tensor 
order parameter. This immediately implies eq. ([L2|), with 



and 



XI = (2/3)Axo, 



Ax^2AxoQo{T,H). 



(58) 



(59) 



Equation ( pSj ) assures that xi is perfectly smooth 
through TjsfT, since Axo is; this follows from the fact 
that |Sp = 1, regardless of temperature. Hence, taking 
the trace over a and (3 in equation (|5^) fixes the coeffi- 
cient of Saf3 at exactly 1/3, independent of temperature 
(the first term is traceless by construction and does not 
affect this result). This coefficient, in turn, leads to the 
exact relation ( pq ) between xi snd Axo, which in turn 
implies analyticity of xi- 

We must still determine the scaling of Qq (T, H) . Since 
Qo involves correlations of two S fields, its renormal- 
ization group eigenvalue is simply (see the appendix) 
—2Ph/i^h- Therefore, the renormalization group trans- 
formation near the d = 3 Heisenberg critical fixed point 
can be used to relate Qo {T, H) to Qo at rescaled temper- 
ature and field as 



Qo(T,i/) = 6-2'3«/''"go(&'/''"t,6^'^ff), 



(60) 



where yh is the renormalization group eigenvalue of H . 
This eigenvalue can be determined by rewriting the mag- 
netic Hamiltonian as 

7^,„ag = Axo HaHpiS^S^ - 1/3|S|2<5„;3) (61) 



2/r = 2 1 - 



- 18n - 88 



n + 8 4(n + 8)2 
1.77 ±0.01, 



(62) 



where in the second, approximate, equality we have set 
d — 3 and n — 3. The quoted error is typical of second 
order in e expressions for d = 3 0{n) models. Since the 
coefficient of the spin tensor perturbation ( |6l| ) is , the 
renormalization group eigenvalue of H is yh — yr/^ ~ 
0.89 ± 0.01. With this eigenvalue in hand, we can now 
obtain the scaling of Qo- Choosing b = \t\~'^'^ in eq. 
(|60|), and using ( p9| ) to relate Qo to Ax, we immediately 
obtain eq.(|l^) for Ax (Item 7), with 



fiix)^Qoit = ±l,H = x). 



(63) 



To obtain the asymptotic form (|T^) of /+, note that deep 
within the disordered phase the susceptibility of Qo to H"^ 
(the coefficient of the magnetic perturbation) should be 
finite. Hence Qo oc in that regime. (Qo of course 
vanishes as 7J ^ in the disordered phase, so this 
term is the leading order term.) In the ordered phase, Qo 
remains non-zero as — > 0, which implies the asymp- 
totic form (|l|) for /_. Finally, requiring that Qo(i, -ff) 
remain finite and non-zero as i — > for fixed H implies 
the final asymptotic form in eq. (^5|) for f±{x) in the 
critical {x —>■ oo) regime. 

All the other results eqs. (|l5-18) quoted for the mag- 
netic susceptibility now follow straightforwardly from the 
scaling law eq. (13) and the asymptotic expression in eq. 
(Ill) for /+. 



VI. FREEDERICKSZ INSTABILITY 

Regarding the Freedericksz instability, we need only 
derive the scaling form (|2^) for Ki, and the asymptotic 
forms (|2^). The former can be done in a manner precisely 
parallel to the argument given for Qq. The difference 
is that the Frank constants have renormalization group 
eigenvalue 1, as can be seen by requiring that the zero- 
field Josephson relation Ki oc \t\'^ holds. 

The first asymptotic form in eq. (|2^) then follows from 
the standard result Hp oc L^^, which holds when ne- 
matic order is well-developed. The second expression in 
eq. ( p2|) follows by requiring that Hp remain finite and 
non-zero in the critical regime. In the disordered phase, 
boundary effects surely fall off exponentially with the dis- 
tance from the boundary for distances large compared to 
the correlation length. This gives the last (exponential) 
form in eq. (p^). 



which differs from (pq ) only by a constant (independent 
of S), owing to |S|^= 1. In this form, Timag is simply 
the "spin tensor interaction" of reference [Q . The renor- 
malization group eigenvalue of this perturbation has the 
e = 4 — d expansion Q 



VII. MAGNETIC SUSCEPTIBILITY NEAR I/T 
TRANSITION 

Finally, we consider the magnetic susceptibility and 
depolarized light scattering near the // T transition. All 
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of the manipulations up to equation ( p9| ) are the same 
here as near the N/ T transition. The differences all arise 
in the scaling behavior of Qo{t, H). This can be obtained 
by noting that, since there is no spontaneous nematic 
ordering anywhere near the // T transition (except at the 
multicritical point), the nematic director will necessarily 
be along the applied field H. (Without loss of generality, 
we will consider H along the z-direction.) In that case. 



(64) 



whence 



Qo = ^[{iS[f)-l/3]. (65) 
Since the magnetic Hamiltonian can be rewritten 



H 



mag 



= Axoi?'V[(Sr)'-l/3] 



by adding a trivial constant, it follows that 

3 dF 



2Nd{h^) 



(66) 



(67) 



where N is the number of sites in the system, and h = 

Using (|67|), we can obtain Qq by familiar manipula- 
tions of the partition function. When Timag is added to 
the original gauge theory Hamiltonian, the resulting par- 
tition function is 



l'Y[dn,J2cxplKj2u{P) + 

C/,,S, . S, + /i2^[(5f )2 - 1/3]|. (68) 

As before, we perform a "polymer" expansion of the 
terms g^iJUijSiSj) -j-j^g partition function. Extracting 
a factor 



(69) 



which represents completely decoupled spins {{■) ^ indi- 
cates normalized expectation with completely uniform 
spin distribution on the unit sphere), we rewrite the par- 
tition function as 



Z = ZdZgY, (n^(^'") 

{C„} \ m / 



(n 



^(js,.s,)n 



(n(ji)ec™ n 



(70) 



This is a sum over disconnected polymers C„, with fac- 
tors of e'^ [^■^•'^ -1/3] decorating the polymers; the same 
factors for sites not on polymers have been absorbed into 
the Zd out front, which requires the denominator to 



avoid double counting. Zd is clearly a smooth function of 
H. Near the // T transition the singular contributions to 
Qo come from the fluctuations of the J/^ 's. For small J, 
it is safe to expand in J. The calculation of the effective 
gauge coupling is similar to that in Paper I (section VII 
D), and involves resumming the leading order contribu- 
tions, i.e., the polymers consisting of single disconnected 
plaquettes. The result is 



K,s = K + J' 



(n(S.-S,-)ne^^'(^-)^-^i) 



(71) 



To calculate the lowest order contributions in H and J 
from eq. (|7^), note that the denominator can be replaced 
by unity. In the numerator, the product of • Sj terms 
only depends upon the relative orientation of the spins 
in the polymer. For the order term, the absolute 
orientation of only one spin matters. Integrating over 
the orientation of this spin first, with relative orientations 
fixed, and then over the relative orientations, the result is 
zero since {{Sf)'^)^ = 1/3. Thus, ||^ the expansion starts 
at order 

The coefficient of ff ^ is easily calculated to lowest order 
in J. Eq. (|7l|) becomes 



Kr,ff = K 



176 

6075^ 



(72) 



The most important consequence of this analysis is that 
the generic singularity of the magnetic susceptibility near 
the IJT transition has the 1 — a/ singularity reported in 
eq. ( psi ) (Item 9). 

All of the critical behavior of Qo can be obtained from 
the field dependence of i^cff via equation (B^, which gives 



Qo 



3 dF aXcff 



2N dK^s d{h^) 



3^ 

' N 



C(J) 



dF 



= (a + 6|t|i-")ij2 



(73) 



where in the last equality, we have inserted the known 
singularity of F at the (zero-field) // T transition, as well 
as the fact that dF/dK^s oc dF/dT. As a bonus, ^) 
also tells us the field dependence of the transition tem- 
perature TjT- 



Tit{H) = Tit{0) + const, x 



(74) 



VIII. DEPOLARIZED SCATTERING NEAR I/T 
TRANSITION 

This hard-won knowledge of the singular behavior of 
Qo near Tjt enables us also to determine the critical 
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behavior of the depolarized hght scattering I^z (q, t) near 
Tit (Item 3). As shown earher, its Fourier transform 
Ixz{r,t) obeys 

/.,(q)cx(Q,,(r)Q,,(0))^a,(r). (75) 

We can now obtain the behavior of (r) from scahng: 

Cy(r,t) = b^''^ajib-^r,b'/''t), (76) 

where Xq is the renormahzation group eigenvalue of Qij 
and b the rescaling factor. We can extract Xq from the 
results just obtained for the singular part of the non- 
linear susceptibility. Scaling implies 



(77) 



where Qq'"® is the singular part of the magnitude of the 
nematic order, and Xh is the (as yet unknown) renormal- 
ization group eigenvalue of H at the Ising gauge critical 
point. Choosing b — \t\~'^, we obtain 



QrHt,h)^\t 



h 



where 



r™(x)^Qf'(l,x). 



(78) 



(79) 



For this to be consistent with our earlier result (|7^) 
for Qo, we must have f{x) oc at small x and also 
|^|-i'(js:Q+2Ah) ^ which implies 



Xr 



2Xh. = 



a-1 



(80) 



To obtain Xq, we still need Xh- But since couples 
directly to Qo {i.e., to S"^ ~ 1/3, whose average is Qo), 
we know from scaling that 2Xh — d + Xq. Making this 
substitution, we find (with d set to 3) 



(81) 



Using this eigenvalue in ([Tq), and choosing b = \t\ " , we 
obtain the scaling form 

a,(r,t) = |tp''+i-"/,,(r/C), (82) 

where ^ oc \t\~'^ is the Ising transition correlation length. 
As the critical point is approached {t — > 0), Cij{r) re- 
mains finite but non-zero at fixed r, so that the function 
fij must scale in such a way as to offset the t dependence 
of the prefactor in This gives the result 



Cy(r,t = 0) ocr("~i'/''-3, (83) 
which, when inserted back into (|7^), implies 



Ixz{q,t = 0) oc 



|.3— (a — l)/z^ 



02 + 62|q 



{l-a)/u 



(84) 



where 02 and 62 are constants and the exponent {l — a)/^ 
has the numerical value 1.44. Note that the scattering 
remains finite at |q| — 0, and lacks even a cusp since 
the exponent q is greater than one; two derivatives are 
needed to get a divergence {d^Ixz/dq^ ~ |q|~°'^^)- 

Finally, we argue that a lack of microscopic head-tail 
symmetry will not alter these results near the // T tran- 
sition. It suffices to show that the leading order contri- 
bution to KcS remains 0{H'^), even in that case. The 
term is non-vanishing (even a system lacking head- 
tail symmetry will still have a Axo^(H. ■ S^)^ term in 
its magnetic energy), and we need show only that lower 
order terms continue to vanish. This is an immediate 
consequence of the absence of nematic order, even in the 
presence of such a microscopic symmetry breaking. No 
analytic scalar (rotationally invariant) function can be 
made which is linear or cubic in H, and the term 



vanishes by equation (67) 



IX. CONCLUSION 

In this paper, we have analyzed critical behavior of 
a variety of quantities at the transitions into the new, 
topologically ordered phase of nematics predicted in our 
gauge model |jl|] . Characteristic signatures were found for 
several standard experimental probes of soft condensed- 
matter, notably specific heat, and optical and magnetic 
response. We hope that these behaviors will be sought 
in the laboratory in the near future. 
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X. APPENDIX 

In this appendix we calculate the scaling dimension of 
Qo{T). Writing the spin as S = M{t)z. + 5S, M being 
the magnetization arbitrarily assigned to point along z. 



(Q,,) = M\t)S,zSjz + {6S,6Sj - %/3). 



(85) 



The second term is the expectation valud of the con- 
nected piece of the operator r = SiSj — {l/3)Sij. This 
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expectation vanishes as t^'^ as t ^ 0. By general renor- 
malization group theory Q , 

= Krf-2/r), (86) 

where yr is the RG eigenvalue of the operator r at the 
relevant (i.e. n = d = 3 Wilson-Fisher) fixed point. 

Thus, the second term in (|8^) is asymptotically domi- 
nated by the first as < ^ provided 2/3 = z/((i + 7^ — 2) < 
/3r, or equivalently, if and only if 

2(3 -(3r = HVr + ?7 - 2) < 0. (87) 

The epsilon expansion for j/^. can be found in Ma's book 
(ref. g, pg. 355), and is 

^-^%^-^+ 4(. + 8)3 -^j+^(-^)' (88) 

which, combined with rj = 2(n+8) ^ ^{^) at n = 3 
gives 

which is indeed negative, as desired. The expectation 
( ^ ) is dominated by the connected piece and scales with 
an exponent 2/3. 



namely mo^.H ■ Si, where mo in a constant perma- 
nent molecular dipole moment, which would dominate 
the quadratic (in H) term (^). For this to make sense, 
we have to rewrite it in a gauge invariant form just as 
in Paper I with the likewise gauge-non-invariant ■ Sj 
perturbation. To wit, we introduce additional variables 
CTi = ±1 on the sites which keep track of whether S; 
points along or opposite the local dipole, and which trans- 
forms like S under a gauge transformation so that the 
combination Oi Si is invariant. Then the perturbation be- 
comes mo^.((TiH) • Si. Since cr's on different sites are 
decoupled from each other, it is trivial to sum over them, 
yielding a magnetic perturbation of 

"Hmag = ^ In [cosh(moH ■ SO] , (90) 

i 

an explicitly even function of each Si. Expanding to low- 
est order, we recover the earlier expression (p5|) with 
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